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Relaxational effects in stellar heat transport can in many cases be significant. Relativistic Fourier- 
Eckart theory is inherently quasi-stationary, and cannot incorporate these effects. The effects are 
naturally accounted for in causal relativistic thermodynamics, which provides an improved approxi- 
mation to kinetic theory. Recent results, based on perturbations of a static star, show that relaxation 
effects can produce a significant increase in the central temperature and temperature gradient for 
a given luminosity. We use a simple stellar model that allows for non-perturbative deviations from 
staticity, and confirms qualitatively the predictions of the perturbative models. 
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1 INTRODUCTION 

The formation of compact astrophysical objects such as neutron stars and black holes is usually preceded by a period 
of radiative gravitational collapse in which neutrinos and photons carry energy to the exterior spacetime. The problem 
of radiative collapse was originally addressed by Oppenheimer and Snyder (1939), who investigated a collapsing dust 
sphere. Their idealised treatment has subsequently been improved via more physically realistic thermodynamical 
models. 

Most work has been based on the standard Fourier-Eckart theory of heat transport (see, e.g., Grammenos, 1994). 
However, this theory is non-hyperbolic, and thus suffers from the pathologies of infinite propagation speed of thermal 
signals, and instability of all equilibrium states (Israel and Stewart, 1979, Hiscock and Lindblom, 1983). The causal 
(hyperbolic) generalization of standard thermodynamics, due to Israel and Stewart (1979), is based on an improved 
kinetic theory approximation, and avoids these pathologies. The causal behaviour results essentially from treating 
the heat flux as an independent dynamical quantity, on the same footing as the equilibrium quantities. As a conse- 
quence, non-quasi-stationary and transient relaxational effects are incorporated, expressed via the appearance of a 
time derivative of heat flux in the transport equation. By contrast, the standard quasi-stationary transport equation 
contains no time derivative; it is an algebraic relation for heat flux, incapable of accounting for relaxational effects. 

In relativistic stellar collapse, strong deviations from quasi-stationary heat transport can occur, so that relaxation 
effects can be significant, as argued by Di Prisco et al. (1996) and Herrera and Santos (1997). Recently, a number 
of radiating collapse models with causal heat transport has been developed (Martinez, 1996, Di Prisco et al., 1996, 
Herrera and Santos, 1997, Herrera et al., 1997, Govender et al., 1998, Herrera and Martinez, 1998).[] 

These models show that the causal theory can predict significantly different results when the relaxation time-scale 
is sufficiently large relative to other time-scales. Typically, for a given luminosity, causal corrections lead to a greater 
temperature gradient, and thus greater central temperature. In other words, non-quasi-stationary heat transport 
implies that a collapsing star with a given observed luminosity is hotter, and with greater temperature gradient, than 
predicted by quasi-stationary heat transport. 

The role played by the relaxation time during radiative gravitational collapse was investigated by Herrera and 
Santos (1997), assuming that the spacetime metric is close to a static (and non-radiating) initial metric. Perturbation 
of the initial configuration is assumed to produce collapse without shearing. In this paper, we use a very simple 
shear-free collapse model with the same static limit as in Herrera and Santos (1997), but which allows for large 
deviations from the initial static metric. Although these deviations are very simple in form, they do allow us to 
provide a non-perturbative complementary approach to the perturbative model of Herrera and Santos. We are able to 



1 See also Triginer and Pavon (1996) and Maartens et al. (1998) for cosmological models with causal heat flux. 
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find explicit analytic forms for the temperature for a simple relaxation time, using the methods of our previous work 
(Govender et al., 1998). We also consider the case of small relaxational effects, to find a perturbative solution for the 
temperature (without assuming that the metric is close to static) . Our results confirm the increase of temperature 
and temperature gradient for a given luminosity. 



2 SIMPLE MODEL 



A shear-free isotropic stellar fluid has metric 

ds 2 = -A 2 (t,r)dt 2 + B 2 (t,r) [dr 2 + r 2 dtt 2 ] , 

in comoving coordinates, where dil 2 is the metric of the unit 2-sphere. The fluid 4-velocity is u a — A~ 1 6 a o, and its 
volume rate of expansion O = V a u a is 

e =!f- <» 

where an overdot denotes d/dt. Since the shear rate a a b vanishes, the expansion rates in the radial and tangential 
directions are the same, and equal to B/(BA). If the star is collapsing, then of course B < 0. 

Modelling the star as a perfect fluid with heat conduction (where the heat flux is provided principally by neutrinos), 
the energy-momentum tensor is 

Tab = (P + P)u a U b + PQab + QaUb + QbUa , (2) 

where p is the energy density, p is the pressure, and q a = qn a is the heat flux, with n a a unit radial vector. 
Herrera and Santos (1997) assumed a shear-free perturbed metric of the form 

ds 2 = - [A {r) + ea(r)Y{t)} 2 dt 2 + [B {r) + eb(r)Y(t)} 2 [dr 2 + r 2 dfl 2 ] , (3) 

where Aq and Bq describe the initial static star, and e <C 1 is the perturbation parameter. It follows that (Bonnor et 
al., 1989) 

p(t,r)=p (r)+eE(t,r), 
p(t,r)=po(r) + eP(t,r), 
q{t,r) = eQ(t,r) , 

where subscript zero denotes equilibrium quantities. The dynamical stability of this particular model has been studied 
extensively in the past by Chan et al. (1989) and Herrera et al. (1989), but a systematic thermodyanamical treatment 
of the heat flux and temperature gradient was not carried out. 

Here we consider an alternative approach to that of Herrera and Santos (1997), where we use a simple metric that 
accomodates large (though restricted) deviations from staticity, and that leads to simple analytic expressions for the 
important quantities. We use the model first proposed by de Oliveira et al. (1985, 1986): 

ds 2 = -A 2 (r)dt 2 + B 2 (r)f 2 {t) [dr 2 + r 2 dn 2 } . (4) 

It follows from equations ([!]) and (Q) that the expansion rate is 

Einstein's field equations for the metric (^) and the energy-momentum tensor (||) give (Bonnor et al., 1989) 
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where a prime denotes d/dr 7 and the pressure isotropy condition 
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In equations (^) and (ffl) we have introduced 
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which define the energy density and pressure of the initial static configuration before the onset of collapse. The 
functions Ao{r) and Bo(r) must satisfy the isotropy condition (JsJ) in order to obtain a solution to the Einstein field 
equations with isotropic pressure. Failure to check this condition can invalidate a solution (see e.g. Govinder et 
al., 1998). Equations (Q)-(|j) together with ( fl0| ) and ( pi] ) describe a radiating shear-free star with an initial static 
configuration. 

The interior metric must be matched to the outgoing Vaidya metric across a time- like hypersurface E, described 
by r = r E = constant in comoving coordinates. The Vaidya (1953) metric is 



ds 2 
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where m(v) is the mass function. The static solution, before the onset of collapse at t = — oo, matches the exterior 
Schwarzschild spacetime across S. In the static limit, the pressure po vanishes at the boundary: (po) s — 0. At late 
times the model becomes nonstatic and the pressure (p) s is nonzero in general because of the presence of heat flow 
(Santos, 1985). The junction conditions resulting from the matching of (^) to the outgoing Vaidya metric yield 
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The proper radius of the star is 
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Using the junction condition ( p. 3D together with (|7|) and (g), and taking into account (p ) 
equation 
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0, we obtain the evolution 
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ispositive because the static solution (Aq, Bo) must match with the exterior Schwarzschild metric. A first integral of 
& is given by 



where the constant of integration is —2ab. We choose b = 1 so that / — ► when / — » 1 in the static limit. Thus 
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Note that p s is nonnegative, so on using the result (po) s = 0, together with equations ([?]) and (|l7|). we obtain / < 0. 
By equations (0) or (fig), this implies that the only possible evolution is contraction. Also, equation ( |l9| ) gives |/| < 1, 
and without loss of generality we assume / is non-negative, so that 



< f(t) < 1 



Integrating equation ( |l9| ) we obtain 
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where the constant of integration has been absorbed by rescaling t — ► i + constant. This shows that / decreases 
monotonically from 1 at t — — oo to zero at t = 0, where a physical singularity is encountered, since /? — > oo by 
equation (^J). From equation dl5| ) we can obtain the simpler expression 
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since at t = — oo the model is static. In the infinite past we must match the static perfect fluid solution to the exterior 
Schwarzschild solution in isotropic coordinates and we obtain 
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where 

is the initial luminosity radius. We can rewrite equation (|l^) as 

a = — ^ , 
fto 2 

where too is given by (p3|). The surface redshift z s is 
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We require q > 0, which implies that 
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For a physically reasonable model the energy density and pressure must be decreasing functions of the radial coordi- 
nate. From equations (Pj and (Q) we can write the spatial gradients of the energy density and pressure as 
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From these equations it follows that in order to have p' < and p' < 0, we must have po' < and po' < for the 
initial static configuration. It has been shown by Bonnor et al. (1989) that the strong energy condition is satisfied 
throughout the collapse if the initial static configuration satisfies the condition 

Po ~ 3po > ^2 ( 31 ) 

at r = 0. 



3 CAUSAL HEAT TRANSPORT 



We now investigate the evolution of the temperature profile using the causal thermodynamics developed by Israel 
and Stewart (1979). The causal relativistic transport equation is 

Th a b u c V c q b +q a = -K {h a b V b T + Tu b V b u a ) , (32) 

where h a b = S a b + u a u b is the projection tensor into the comoving rest space, T is the local equilibrium temperature, 
K is the thermal conductivity and r is the relaxational time-scale. Setting t = in ( |32"| ) we obtain the non-causal 
Fourier-Eckart law for quasi-stationary heat transport, which leads to pathological behaviour in the propagation 
velocity of thermal signals. 

For the metrics (||) and (0), equation (|2|) reduces to 

Tq + Aq = ~{AT)' . (33) 

Taking r ~ 10~ 4 s, which corresponds to the early stages of the evolution of a neutron star, Herrera and Santos 
(1997) found that the magnitude of the temperature gradient was enhanced by 10 per cent as compared to the case 
t = 0. For a temperature of approximately 10 6 K and a thermal propagation velocity v rj 10 3 cm/s, Di Prisco 
ct al. (1996) obtained r i=s 10 2 s. Using this value as an upper bound on the relaxation time, Herrera and Santos 
found that the temperature gradient is five orders greater than the corresponding non-causal case. In general, they 
concluded that for a given luminosity, the associated temperature gradient is enhanced by relaxational effects, and 
the enhancement grows with increasing relaxation times. (This qualitative conclusion is confirmed by Di Prisco et 
al, 1996, and Govender et al., 1998.) 

In this section we perform a similar analysis but with the simple non-perturbative metric (0). For a physically 
reasonable model, we assume that heat is carried to the exterior spacetime by thermally generated neutrinos having 
long mean free paths through the stellar core. The thermal conductivity is then given by the radiative transfer form 
(Martinez, 1996) 

« = lT\ c , (34) 

where 7 (> 0) is a constant and r c is the mean collision time for neutrino-matter interactions. These interactions are 
primarily due to electron-neutrino scattering and nucleon absorption. Following Govender et al. (1998), we assume a 
mean collision time of form 

rc=(f)r-', (35) 

where a (> 0) and a (> 0) are constants. For thermal neutrino transport, Martinez (1996) argues that a w §. The 
mean collision time decreases with growing temperature except for the special case a = 0, when it is constant. The 
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special case of constant r c was investigated by Govender et al. (1998) for acceleration-free gravitational collapse. By 
adopting a similar approach, we investigate the relaxation effects when the acceleration is nonzero. The relaxation 
time-scale in causal radiative transfer (Udey and Israel, 1982, Schweizer, 1988) is taken as equal to the mean collision 
time. A simple generalization 



a 



(36) 



where (3 (> 0) is a constant, allows both for cases where the relaxation time-scale is greater than the collision time- 
scale (Pi /a > 1) (Maartens and Triginer, 1998), and for the case of perturbative deviations from the quasi-stationary 
case, i.e. (3y/a <C 1. The degenerate case (3 = recovers the non-causal transport law. 
With these assumptions, the transport equation (33) implies 
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where f(t) is implicitly determined by equation (|2l|). Equation ( |37| ) is a first order radial equation in AqT, which 
is solved in principle once An(r) is specified and boundary conditions are given to determine the time-function of 
integration. The effective surface temperature as determined by an external observer is given in general by Misner 
(1969); for the metric (|) this leads to 
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where £ is a constant, and is the total luminosity for an observer at infinity. This can be found using the general 
form in Bonnor et al. (1989): 
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where the surface redshift z s is given by equation (|26|). 

Equations (|3^ ) and ([39|) determine the time- function of integration of the temperature equation (|37|). In the 
non-causal case (3 = 0, the temperature equation is readily solved. Writing T for the non-causal temperature, we find 
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where T(t) is determined from the surface conditions via equations (|26|), 
we find 
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) and (|39"1). For example, if a = 0, then 
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The temperature equation (37) is difficult to solve in the general causal case (3 > 0. However, we can find a 
perturbative solution in the case of small relaxational effects, and we can also find an exact solution in the case a = 0, 
i.e. when the mean collision time is constant. This assumption may be reasonable over short periods of time. 



3.1 Small relaxational effects 



For a qualitative estimate of how relaxation effects change the predictions of the quasi-stationary theory, we consider 
the case of small relaxation parameter (3j/a, and solve equation ( |37| ) perturbatively. We write 
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where the quasi-stationary temperature T is given by equations (j40|)-(f42|). In order to make sensible comparisons, we 
require the luminosities, and thus the surface temperatures, to coincide, so that 

5 E =0. 

To first order, the temperature equation fii7\) gives 
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This linear equation in AqS has solution 
where T is given in equations (40)-(E3). Now it follows from equations (p^|)— (p^|) that 
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which is always negative, since < / < 1 and / < 0. Thus we see that 

S > 0, 

so that small relaxational effects lead to an increase in the temperature into the interior, consistent with the results 
of Herrera and Santos (1997). 

For the case of thermal neutrino transport, where a = |, we can give explicit analytic forms for the relaxational 
temperature correction. By equation (Ei 
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We find T from the surface conditions via equations (|26|), (|38| ) and (|39| 
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where is given by equation ( [45| ) and F(a, 6, c, z) is the hypergeometric function. Thus, collecting the above results, 
we have the explicit analytic form for the causal temperature correction in the case of thermal neutrino transport 
with small relaxation timescale: 
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3.2 Constant collision time approximation 



For a = 0, equation (|37|) integrates to give 
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where T(t) is a function of integration. Using equations (|46| ) 
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Hence we can finally write the temperature as 
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If we set /? = then (|47j) becomes 
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where T is the noncausal temperature. This is exactly the temperature obtained by de Oliveira and Santos (1985) 
when considering shear-free stellar models in the noncausal theory. Our expression for the temperature in (|47|) differs 
by the additional term involving f3. 



4 CONCLUSIONS 



Using a very simple shearfree stellar metric, we have investigated relaxational effects on the temperature in a 
radiating collapsing star in the non-perturbative regime relative to the static limit. Our results are in this sense a 
non-perturbative (though restrictive, due to the simplicity of the metric) complement to the perturbative results of 
Herrera and Santos (1997), whose model necessarily remains close to the static limit. 

In the case of small relaxation time-scale but possibly large deviations from staticity, we were able to confirm in 
section 3.1 the increase in interior temperature that results from causal heat transport. In other words, heat transport 
that respects relativistic causality requires that the interior of a star must be hotter than predicted by quasi-stationary, 
non-causal transport, for a given luminosity. We also found the explicit analytic form of the temperature correction 
in terms of hypergeometric functions. 

In the case of unrestricted relaxation timescale (section 3.2), we were able to find analytic forms only for the 
constant collision-time regime. A detailed analysis of the causal temperature T in (E7I) is not easy because of the 
presence of the arbitrary functions Aq and Bo. However it is possible to make a few general qualitative statements 
about the behaviour of T. It is easy to see that the causal and noncausal temperatures are different in the interior of 
the star. 

In particular we note that the causal temperature is higher for a given luminosity than the noncausal temperature 
for each interior point of the star since 

4/ -9V? + 5 > 0. 
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Here we have used the fact that < / < 1 and the result 



V 

which follows since Aq > in the interior of the star. The parameter (3 determines the deviation from the noncausal 
temperature T. For small values of /?7 /a, the causal temperature is similar to the noncausal temperature T, but as 
relaxational effects grow for larger values of f3j/a, T and T deviate substantially. 
The causal and noncausal temperatures coincide on the surface of the radiating star: 

COa = (T) s ■ 

By equation ( ff7|) we note that the surface temperature T s becomes zero at t = — oo and at the time of the formation 
of the horizon given by 

/ 2m 
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This is also the case for the noncausal theory, from equation Q471 ) . Note that 

k(T)(A TY - K(f)(A fy = -(Bofr)q . (49) 



Using equations (34) and ( p5| ) we can rewrite equation (J4S|) as 

4-cr 



(T 4 -)' - (T 4 -J = -—^ (A ^B fr)q, (50) 

where we have defined T = AqT and T = AqT. We can immediately see that the relative spatial gradient of T is 
greater than that of T, since / and q are positive. Equation ( |50| ) is similar to the result of Govender et al. (1998). 
However in that case the acceleration was vanishing and the model had a Friedmann limit. In the present model there 
is nonzero acceleration and the model has a static limit at t = — oo. 

To conclude we observe that the line element (g) used by Herrera and Santos is more general than (||). However, 
we were in a position to investigate the evolution of the temperature profile analytically. Using the causal transport 
equation for heat flow we were able to demonstrate the role played by the relaxation time in radiative gravitational 
collapse, providing a non-perturbative support for the results of Herrera and Santos. Looking to the future, a more 
realistic analytic model of dissipative gravitational collapse would need to incorporate the effects of viscosity and 
shear. It would be interesting to investigate the relaxational effects of such a model especially during the late stages 
of collapse when the temperatures can be of the order of 10 9 K (Shapiro and Teukolsky, 1983). 
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